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^ I Abstract 

For a large class of systems of o.d.e.'s which have first integrals, 
^^ ' the method of arrays yields the following results: 

f~^ I i) The first integrals / can be found by solving systems of linear 

equations. 

ii) How the first integral / responds to changes in the system S. 

iii) An easy way for finding the first integral for a special class of 
first integrals if they exist. 



X. 

S 1 Introduction. 



Let S be an autonomous system of o.d.e.'s of the form: 

S: y^ = y^j2c^Jyl'"■■■yn'" (^=l,...,n) 1.1 



In the study of the problem: how a trajectory determined by S, responds 
to changes in the system S; it becomes apparent, that limiting our attention 
to trajectories that are given by first integrals of the form: 
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^efc(y^...y^), 1.2 

fc=i 



the problem can be simplified. For, we can focus on the effect that changes 
in the exponent vectors Hj = {hji, . . . ,hjn) and or the coefficient vectors 
Cj = (cij, . . . , c„.)*, (j = 1, . . . , r), of the system 5", have on the coefficients 
Cfc and the exponent vectors Bk = {bki, ■ ■ ■ , hkn) {k = 1, . . . , g) of the first 
integral /. 

This was done in my paper: Integrals of Multinomial Systems of Ordi- 
nary Differential Equations (Journal of Pure and Applied Algebra, 45 (1987) 
225-240). 

To facilitate the computation, a more efficient notation was employed. In 
this notation the system S of ( |1 . 1| ) becomes: 

r 

S: y' = yY, CjY^^ , where 1.3 

i=i 

y is the column vector (yi, . . . , y„)*, 

Cj is the column vector (cy, . . . , Cn^Y and 

Hj is the row vector {hji, . . . , hjn)- 

Similarly, (|L^) becomes: 

I = J2ekY^\ 1.4 

fc=i 

where Bk is the row vector (6^1, . . . , hkn)- 

We will refer to ( p..3|) and ( p.. 41) as the multinomial vector form of S and / 

(m.v.f. hereafter). 

In this notation, the formula for the derivative of a monomial Y^ , along 
the trajectory, is easily shown to be: 

{Y^y = Y''J2{B-C,)Y''\ 1.5 



where {B; Cj) = J2i=i ^i'^ij (^^e inner product of the vectors B and Cj). Using 
( |1.5| ), the derivative of / in ( |1.4| ) becomes: 



k=i j=i 

neglecting all terms for which [Bk] Cj) = and grouping the coefficients of 
equal monomials together, we get: 

^' = E( E ek{Bk;C^))Y^^ = 0, 1.6 

i = l Ha,+Bk=Ei 



from which the following was proven: 

a) Each B^ (k = l,...,q) satisfies a system of r linear 
equations. 

b) The difference of any two exponent vectors of / is a 

linear combination of exponent vectors of S. 1.7 

c) ei,...,eq is a solution of p linear homogeneous equa- 
tions (the equation are the coefficients of Y^') which must 
vanish, since yi, . . . ,yn are independent. 



To make it possible to classify first integrals by the various relations 
between S and / implied by ( |1 . 7|) , the 'method of arrays' was developed. 
An integral array corresponding to a given first integral / of a system S, is 
a pictorial representation of all the conditions required for J to be a first 
integral of S. 

In this paper we derive similar results for the cases where / is given by: 

J = ln(y^i) + ^efcy^'= 1.8 

fc=2 



/ = eiF^^ + In 1 + ^ efcF^M , 1.9 

V k=2 ) 



with some surprising results. In the case of (|1.8|) , we find that the logarith- 
mic integral is closely linked to an algebraic integral of the form ( |1.4| ) , in the 



following sense. Let /, given by ( |1.4| ), be a first integral of a system 5" given 
by ( |1.3| ), where the coefficients and exponents of S are real or complex num- 
bers. We will show the existence of a multinomial system S{Q\ depending 
on a set of parameters 6 = 6i, . . . ,6m satisfying the following: 



a) S{9) has a first integral I{0) such that euiO) and the 
components oi Bk{6) {k = 1, . . . ,q) are rational functions of 

b) There exist ^such that S{e) = S, 10) = I 

c) If for some value 6*, I{6*) reduces to a non-zero constant, 1.10 
then there exists /*(r) = In (y^i(^*)) + EL2 Cfc^'^"^^*^ 

which is a first integral of S{6*). Thus S{6) which has an 
algebraic first integral given by ( |1.4| ) and S{9*) which has a 
logarithmic first integral given by (|1.8|) , both belong to the 
continuous system S{9) 



In the case where the first integral of S is of the form (|1.9| ) there exists S{6, p) 
where is a set of continuous parameters as above, while p takes on positive 
integral values only. S{9, p) has a first integral given by 

I{9, p) = ei{9, p)Y''^<-''^^ + In ("l + f^ 6^(9, p)Y''>'^''A 1.11 

V k=2 ) 



where 6^(6', p), -8^(6*, p) are rational functions oi 9i, . . . ,9m, p and there exist 
9, p = q s.t. S{9, p) = S, I{9, p) = I. Note that p is the number of monomials 
in/(^,p). 

In addition, the method of arrays yields some curious results, such as: 



i) Let 5* be given by 



n-l 



where Y^^ = J] (y^'^) '\ J = 1, • • • , s. 

i=0 

Let S have a first integral 

g n-l 



/ = E e.i^^M^^^ = n (?/^^0 



then, the exponents of /, Bxik = 1, . . . , g) are independent of lj{j 
1, . . . , s) provided: 



yAf, _^|^^W^)-i^(myn-i) 



1.12 



e.g. Let 5* be given by 

y"=-l,y-'{y'f + l,y + l,y^ 

F*^i violates ([TI^ ), while y*^2,yA^3 (^q ^ot. 5" has a first integral 



/i, the coefficient of y*^^ appears in every exponent of /, while /2,/3 
appear in none. 

ii) Let S he given by 



S : y' = y"^ CjY \ let 5 have a first integral which is either: 



1 
a) I =J2 ^kY^" or 

A:=l 



-Bk 



fc=2 

Let cFa^S) denote the system: 

r 

y' = y"^ CjY" ^ , a any real or complex number. 
Then (Tq(S') has a first integral: 

k=l 

where S^. = aB^, Ck = Ck {k = 1, . . . ,q). 

2 The Derivative Formulas. 



Throughout this paper we assume that 2/1, • • • ,1/n are algebraically indepen- 
dent. Let y = (2/1, ... , yn) be a solution of the system 5* given by ( |1.1|) and 
lety« = yj^..2/^, 

n r n 

{Y^y = Y: hy-'Y^'yl = Y^Y: E(^^c.,)F^^ yielding 



i) of ( |2.1| ) has two interesting consequences 



2.2 



a) A necessary and sufficient condition for a monomial Y^ to be 

a ffist integral of S is that {B; Cj) = (j = 1, . . . , n), this 
implies: 

i) monomial ffist integrals of S are independent of the 
exponent vectors Hj oi S (j = 1, . . . , r). 

ii) Let the matrix (cij) of S have rank s < n, then S has 
exactly n — s independent ffist integrals. For the system of 
homogeneous linear equations: 
(i?;C,) = (j = l,...,r) 
has exactly n — s independent solutions. 

b) Let r = n and let (i/j, Cj) = for all i ^ j, then 

{Y^^y = iH,,Q)Y'^^ (^ = l,...,n). 

Setting Zi = Y^"- yields a complete separation of variables, in 

^1 ) • • • ) Zfi . 



Remark: The sufficiency condition of a) applies to any system of the form: 

r 

y'i = yi^c^jfj{yi,...,yn)- 



b) applies to any system of the form: 



H, 



3 The Method of Arrays. 

Before we give a systematic treatment of the "method of arrays" , we give a 
few examples to show how to construct the continuous system S{9) and its 
first integral I {9), to which a given system S and its ffist integral / belong. 



This can be done for any system S of the form (LI) whose ffist integral 
/ is such that 

r = ^^ where A{y),B{y) 

B{y) 



are linear combinations of monomials. 

Example 1. 
iy'i = y2 = yi{yi'y2) or y' = y{C^Y"^ + C,Y"^) 

\ 1/2 = -yi = y2{-yiy2^) 

where, 

H, = {-1,1), //2 = (1,-1), C,= (lYc2=(\ 
(the latter is the m.v.f. oi S). S has a first integral: 

I = yl + yl = eiY^' + e2Y^\ 

where Bi = (0,2), B2 = (2,0), ei = 63 = 1 



3.1 



3.2 



Applying the derivative formula i) of ( |2.1| ) to /, we get: 

/' = ei{Yf)' + e2{Y''^y 

= ei[(5i;Ci)y^^+^i + (5i;C2)r^^+^^] 3.3 

+ e2[iB,; Ci)y^2+^i + (S2; C2)r^^+^^] = 0. 



Substituting for Bi, B2,Ci,C2, Hi, H2 their values, given by (p.l|), (p.2|). 



in (|3.3|) we find: 

(5i;Ci) = (S2,C2) = 

Bi + H2 = B2 + Hi = El, so that 

/' = [ei(Si;C2) + 62(^2; Ci)]y^^=0 

Thus, the relations between / and S are: 

2) {Bi;Ci) = {B2;C2) = 

li) Bi + H2 = B2 + Hi = El 3.4 

III) ei{Bi;C2) + e2{B2;Ci)=0 



These relations are summarized by the array: 

Which obeys the following: 

i) Ha appears in the k column of a p x g array if and only 
if (Sfc,C„)^0. 

ii) Ha appears in the j row and the k column oi a p x q 
array if and only if: Bj. + H^ = Ej, 1 < j < p, ^ < k < q. 

In our case the array tells us that: 
/'=[ei(Si;C2) + 62(^2; Ci)]y^^=0 
where 
El = Bi + H2 = B2 + Hi. 

Now, ii) of (^.41) and i) of (|3.5|) imply: 



[Bi-C^) = iBi-B2;C2) + iB2;C2) 

= {Hi - H2; C2) + = {Hi - H2; C2) + 

[B2\Ci) = {B2-Bi-Ci) + {Bi-Ci) 

= {H2-Hi;Ci) + = {H2-Hi;Ci)j^0 



3.5 



3.6 



We can, now, use i) of ( |3.4| ) and ( |3.6[ ) to solve for Si, 52,61,62- For 



Bi, B2 are solutions of the linear systems: 

r {Bi;Ci) = r {B2;C2)=0 

'> \ {Bi- C2 = {Hi - H2; C2) ''' \ {B2; Ci) = {H2 - Hi- Ci). -^-^ 



These systems have unique non-zero solutions provided: 
i) d = det f ^11 ''0^0 

V C21 C22 ; 

ii) {Hi-H2;Ci)^0 
III) {Hi-H2;C2)j^0. 



We can now use in) of ( |3.4| ) to find 61,62- 

Since the only conditions on 5" are the inequahties (|3.8|) , we may take for 
S{9), the continuous system to which S, given by ( |3.1| ), belongs, the full 8 
parameter system: 



S{9): y' = y[ 



Cll 

C21 



-y^(/lll,/ll2) _|_ 



C12 
C22 



y(/l21,fi.22)l 



3.9 



subject only to the inequalities ( |3.8|) . 

We can, now, solve the linear systems (137 



i) B,{e) -- 
It) B2ie) 



^^^^-^(-C2i,cn) 

= ^^^^-f^(c22,-Cl2) 



for Bi{9), 32(6), and obtain: 



3.10 



To check that ( p.lO|) are solutions to (pl7|), note that (— C21, cn), (C22, — C12) 
are normal to Ci, C2 respectively and that 

((-C21, Cll); C2) = d ((c22, -C12); Ci) = d. 

B2{6) may also be obtained, when Bi{6) is known, by using ii) of 
yielding: 



B2{e) = Bi{e) + H2 - Hi 



3.11 



To show that (|3.11| ) agrees with ii) of ( |3.10|) , we show that ( p.ll| ) is a 
solution of the linear system: 

{B2{e);C2)=0 

{B2{e);C,) = {H2-Hi;Ci) 

For 

{B2{e)-C2) = {B^{e)-C2) + {H2-Hr-C2) 

= {H,-H2;C2) + {H2-Hi-C2)=0 
(i?2W,Ci) = {B^{e),C^) + {H2-H^,C^) = {H2~H^,C,), 



10 



Thus ( 3.11 ) and ii) of ( |3.10 ) are both solutions of the hnear system ii) 
of ( p.7|) which has a unique solution when the determinant d is not equal to 
zero. 

To find ei{9), 62 (^) we use Hi) of 



ei(//i - H2; C2) + e2{H2 - H^,Ci) = 0, yielding 

ei = {H2-Hi;Ci), 62 = {H2 - Hi,C2) 

Thus 



3.12 



C22„,-Cl2 



+ {H2-H,,C2)iyry2 



(H2-gl.gl) 



setting C12 = C21 = 0, Cn = 1, C22 = -1 we get d 
(-1, 1), H2 = (1, -1) in dp) and (gJ^) we get: 

S{e) = S of Q and I{9) = 21 of Q 



-1 and Hi 



We now show that the system given by ( p. 9]) has a first integral even when 
the inequalities are violated (one at a time). 

Let d = then C2 = ICi and by i) of (2.2) 5* has the monomial first 
integral Y^ where {B; Ci) = (5; C2) = 0. Now, let rf ^ and let 6* be such 
that {Hi - H2, C2) = while {Hi - H2, Ci) 7^ 0, then 

I{9*) = {H2 — Hi, Ci) and fails to define a first integral of S{9*). Fortu- 
nately, {Hi — H2, C2) = is the very condition required for S{6*) to have a 
first integral of the form: 



rid*) = in y^^ 



yB2 



3.13 



For applying the derivative formula ii) of (^TTj) to I* {9*) we get 

{r{9*)y = {Bi; C2)y^^ + {B2; Ci)y^^+^^ = 3.14 

Thus the relations between S{9*) and I* {9*) are: 
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i) {B,;Ci) = {B2;C2) = 

ii) H2 = B2 + Hi = El 3.15 

m) (5i;C2) + (S2;Ci) = 

ii) of ( |3.15| ) implies B2 = H2 — Hi and i) and Hi) of ( |3.15| ) yield the system 
of linear equations: 



[Bi,Ci) = 

[Bi, C2) = — (-B2, Ci) = {Hi — H2; Ci] 



A n {Hi- H2\Ci) 
and Bi = (-C2i,cii) 



3.16 



^' f''^'^ In {vi'^'^yr) + yf--'^"?/^-'^- 3.17 



The integral array of /*(6'*) is A* = {H2 Hi) (the o above H2 indicates that 
Bi is not added to H2 to get ii) of ( ^.151 )). Similarly, if S{6**) is such that 
{Hi — H2, C2) 7^ 0, d 7^ but, {Hi — H2] Ci) = then the relations between 
r{9**) and S{e**) are: 

i) {Bi,Ci) = {B2,C2) = 

ii) Bi + H2 = Hi = El 3.18 

m) {Bi;C2) + {B2;Ci) = 

the integral array of I* {9**) is 
A** = {H2 Hi) and 



^(.n-/..)^(/^i.-/^..) + {H2-H1.C2) j^ (^c..^-c,.) 3_^g 
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Summarizing the above: 

The system S of ( |3.1|) and its first integral / of (|3.2| ) belong to the con- 
tinuous 8 parameter system S{9) of ( |3.9[ ) and its first integral I{9) of ( |3.12| ) 
which exists provided S{9) satisfies the 3 inequalities: 

i) d^O 

ii) {H,-H2;C,)y^0 

in) {H^-H^-C^)^^ 

If d = 0, S{6) has the monomial integral / = Y^ , where {B\ Ci) = {B; C2) = 
0. 

If S{e*) is such that d ^ 0, {Hi - H2, Ci) ^ but {H^ - H2; C2) = 0, 
then S{6*) has the logarithmic integral I* (9*) given by (|3.17|) . 



If d ^ 0, {Hi - H2; C2) ^ but {Hi - H2, Ci) = 0, when 9 = 9**. Then 
S{9**) has the logarithmic integral I* {9**) given by ( |3.19|) . 



4 

Let 5*3 (6') be the system given by: 

y' = y[CiY^' + CaF^^ + CsY^^], where 

Ci, Hi,C2, H2 are as in ( p.9|) , subject only to the three inequalities (|3.8| ). We 
shall refer to the system S{9) of ( p.9| ) as 5*2 (6') and write 

S3{9) = S2{9)+CsY''' 

Since n = 2 and Ci, C2 are linearly independent, we may write 

C3 = liCi + I2C2 4.1 



Now, 5*2 (6') has a first integral l2{9) given by ( p.l2| ). We are going to show 
that 5*3 (^) has a first integral 

/3(0)=/2(^)+e3F^3 

subject only to the following: 
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(ifs; C3) = h{H,; Ci) + /2(^2; C2) 4.2 



where /i,/2 are as in ([4.1|), and Ss{6) satisfies the following additional in- 
equalities: 

r t) iH,-Hs;C\,Cs)^0 



There are, now, 2 cases to consider. 
Case 1. 

hxk^O 4.4 



This inequality and the inequality d = di2 7^ imply that any two of 
Ci, C2, C3 are linearly independent. To find h{6), let 

Bs = B, + H3- Hi 4.5 



where i?i, is as in i) of ( p. 101 ). We show that ( |4.2| ) implies (-B3; C3) = 0. For, 
using; i) of f|3.7|) we get: 
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(53; ^3) = {Bi + H,-Hi;Cs) 

= {Bi;Cs) + {Hs,Cs)-{Hi,Cs) 

= {Bi;liCi + kC2) + li{Hi,Ci) 

+ l2{H2;C2)~h{Hi,Ci)-l2{Hi,C2) 

= l2{Bi,C2) + l2{H2-Hi;C2) 

= l2[{Hi-H2,C2) + {H2-Hi,C2)]-- 





Also 




a) 


(53,Ci) 


= {Bi + H,-Hi-Ci) 

= {Bi;Ci) + {H,-Hi;Ci) 

= + {H,-Hi;Ci)y^O 


in) 


{B3;C2) 


= {Bi + H,-Hi;C2) 

= {Bi;C2) + {H,~Hi;C2) 

= {Hi-H2-,C2) + {H,-Hi-C2) 

= {H,-H2;C2)j^0 


iv) 


(BuCs) 


= iBs + Hi-H,;Cs) 

= iB,,C,) + {Hi-H,;Cs) = {Hi 


v) 


B2 + H, 


= -Bi + H2 — Hi + H^ 

= Bi -\- H^ — Hi + H2 = -B3 + H2 


vi) 


iB2;C,) 


= {B2-B,;C,) + iBs;C,) 






= {H2-H,,C3)+0y^0 



4.6 



^3; Ca) ^ 



To get the results of (pj) we used (|J), (U), (gj), (gj) and 
The relations ( [4.6| ) are summarized by the array: 

Bi B2 -D3 

j^^ , H2 Hi 0\ El 

^ Hs Hi \ E2 ^-^ 

H, H2 Es 



Columns k says [Bk, Cj) 7^ if and only if A; 7^ j {k = 1, 2, 3) Row 1 gives 

the relations (|3.4|) . 

Row 2 implies Bi + H^ = B3 + Hi = E2. 

Row 3 implies B2 + H^ = B^ + H2 = E^. 

Thus, 

/' = [ei(5i,C2) + 62(^2; Ci)]F^^ 
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+ [ei(Si;C3)+ 63(^3; Ci)]F^^ 

+ [62(^2; C3) + 63(53, C2)]F^^ = 



and {Bk, Cj) are given by (|4^ ) (j, A; = 1, 2, 3). 

To find h{(^) of 5*3 (^), we liave l2{0) is as in ( [3.12 ) and we solve for B^ by 

using tlie linear system: 



(B-i, Ci) = {H3 - Hi, Ci] 



4.8 



Yielding: B^ in the same form as Bi, B2 e.g.: 



{H3- Hi,Ci) 
B3 = (C23, -C13), ai3 



d 
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Cll Ci3 
C2I C23 



ci, 62 are as in I2 given by (|3.12 



Setting the coefficient of Y^^ to zero we get 



63 



-ei 



(i^i;C3 
(53; Ci 



-ei 



[H^-Hi-Ci) 



4.9 



We now must show that 

;Si;C2) (S2;Ci) 

;Si;C3) (53; Ci) 

(52; C3) (53; C2) 



4.10 



Now, 

(i?l, C3) = /l(i?l, Ci) + /2(Sl, C2) = /2(Sl, C2) 
(52,C3) = /l(52,Ci) +/2(52,C2) = /l(S2,Ci) 

Thus ( [4.10|) becomes: 



(i?l,C2)(i?2,Ci 



1 1 

k (53; Ci) 
/i (53; C2) 



4.11 
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{B,,C2)iB,,C,)[~hiB,,C,) - l2{B,,C2) 



{B,,C2){B2,C,)[~{B3;C,)] =0 by ^ of H). 
Setting ei = 1 (in hid)), we get: 

e. = ^|^-^#^, es = ^l^^l^i^ by O 



i^i — H2; Ci] 
Thus: 



{H,-H;-C,) 



W) = W) + 



{Hi - H^^Cy, 

{Hi — H^.Ci 



iy?'y2 



(H:i-H^;Ci) 



where l2{0) is the first integral given by ( |3.12| ), multiphed by (jj- _h -c ) ■ 
Let S2{0),l2{9 be as in (|3.1|), (|3.2| ) respectively. Then 



C, 



h 



and 



(^3;G3 



-li + Z2 which implies 



/l(/l31 + l)=/2(/i32 + i: 



4.12 



and B3 = B1 + H3- Hi = {h^i + 1, /iga + 1) subject to (|]T2]). 
{Hs-Hi;Cs) 2/2 



es 



{Hs-Hi;Ci) /131 + 1 



Thus 

/3(^1 = /2(^1 + 63^^^ = yl + yl + ^^^-'^-+i./.3.+i 



hzi + l 



yi 2/2^ 



4.13 



is a first integral of 

^3W = 52(^)+C3y^^ 
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-yi - I2y1''y2''''' 



y'2 



4.14 



subject to ( ]I7T^ ). 

Case 2. One of /i, I2 is zero, say /i, then: 

C3 = ^2^*2 

Let (i^3, C3) be subject to the condition: 
{H^, C3) = /2(^2, C'2) which imphes: 



[H^-H2,C2) = ^ 



4.15 



In addition, let Sz{6) satisfy the following inequalities: 



4.16 



Let Si, B2 be as in ( ^71^ then: 
(5i,C2)7^0 implies (5i,C3)^0 

(^2, C2) = implies (52, C^) = 

Let B3 = Bi + H3- Hi then 

(S3;Ci) = {Bi + H3-Hi;Ci) 

= (5i; Ci) + (if3 - Hi; Ci) = + (//3 - i/i, Ci) 7^ 

(S3,C2) = {Bi + Hs-Hi,C2) 

= {Bi,C2) + {H3-Hi;C2) 

= {Hi-H2;C2) + {Hs-Hi;C2) 

= {H,-H2;C2) = 



by ( [4. 151 ) and 

(53;C3) = (S3,C2) = 

Thus 



= 



All the relation between Bi, B2, B3, Ci, C2, C3, Hi, H2, H3 are summarized by 
the 2x3 array: 



Bi B2 


S3 




112 Hi 


\ 


El 


H3 


Hi 1 


E2 



Bi, B2 are as in hiO) and B^ = Bi + H^, — Hi may, also, be obtained by 
solving the linear system: 

(53,C3)=0 = /2(53,C2) 

(53, Ci) = {H, - Hi- Ci) which yields 

{H^- Hi]Ci) 

B3 = (C22 - C12) 

"12 

The vanishing of the coefficient of Y^'^ yields: 

{Bi,C3) _ -eMBi-C2) 
""' ""'(B.^Ci) (Bs^Ci) 

setting ei = 1 we get 

{Hi-H2,C2) l2{H2-Hi,C2) 



63 = -k 



{H3 — Hi, Ci) {H3 — Hi;Ci 



Thus, if S^{e) = S2{6) + CsY"^, where 32(6) is as in Q, C3 = ^2^2 and 
H3 is subject to the condition: 

(iJ3-^2,C2)=0 

and the inequalities {H3 — Hi; Ci, C3) 7^ 0, then Ss{0) has the ffist integral: 

[H3 - ni,L.i) 
where l2{0) is given by ( p.l2| ). 
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Example 3. 

y" - 2|/'2 + 3y2 = 5.1 



has a first integral 

/ = -% + In \\ - —y''' + 82/2 + 4?/ 

Let y = yi, y' = y2, then the multinomial vector form for ( ^.1|) is: 
2/' = y{CiY^' + Cs^^^ + CsY^^), where: 

ifi = (-l,l), //2 = (0,1), i/3 = (2,-l) 
^^ = ( j ' ^2 = ( 2 j ' ^3 = (^ _3 j • 

J = eiY^' + In (l + 62^^' + esF^^ + e^Y^*) , where: 
5i = (l,0), S2 = (0,2), S3 = (2,0), S4 = (l,0) 

ei = -4, 62 = r-, 63 = 8, 64=4. 



Now, 



.B. V ^ e2(i/^^)- + e3(y^-)- + e4(y^^)- 



'^ ^1 + 62^^2 + 63^^3 + 64(^^4) 
Clearing of fractions and using the relations: 
(52;Ci)=0, {B,,B,,B^;Cj) = ifj^l, 
yields: 
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+ e2[(S2; C2)F^^+^^ + (S2; C3)F^^+^^] 
+ 63(53; Ci)F^3+/^i + e^(5^. Ci)y^^+^^ 







5.2 



Grouping all the coefficients of the same vector together, we get: 

/' = [ei(5i;Ci) + 64(^4; Ci)]F^^ 
+ [62(^2; C2) + eie2(Si;Ci)]F^2 
+ [62(^2; C3) + 6163(^1; Ci)]y^3 

+ [63(53, Ci) + 6164(^1; C,)]Y^^ = 0, 



5.3 



where: 

i) Ei = Bi + Hi = Bi + Hi = (0, 1) 
n) E2 = B-2 + H2 = Bi + B2 + Hi = (0, 3) 
nz) ^3 = 52 + ^3 = 5i + 53 + 5i = (2, 1) 
iv) 54 = 53 + 5i = 5i + 54 + 5i = (l,l) 



5.4 



The integral array A of I is: 



A 



5i 


52 


53 


54 


B1 + B2 


5 


1 + 53 


5i 


+ 54 




( H, 








5i 












\ 


El 





52 








Hi 












52 





53 













5i 







53 


V 





Hi 















5i ^ 


54 



The bars over 52, 53, 54 in columns 5,6,7 are used to indicate that they 



do not appear in the coefficients of Y 



M(A) 



( {Bi-Ci) 

(52; C2) 

(52; C3) 

V (53; Ci: 



= 2,3,4). 

[54; Ci) 








(5i,Ci; 





i) of ( |5.4| ) implies 5i 

ii) of (|5.4|) implies 5i 

Hi) and if) imply 52 

53 



54 = (1,0) 

54 = 52-5i = (l,0) 
352-25i-53 = (0,2) 
25i = 2(52-50 = (2,0) 






[Bi-Ci] 







(5i;Ci) ; 



5.5 
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Setting the coefficients of Y^ = 0, in (|0|), (z = 1, . . . , 4) we get: 

(-Bi;Ci) o (Bi;Ci) -16 



5.6 



We now use the 4x7 array A to find the conditions that <S'(0) has to 
satisfy, to have a first integral 

/(e) = eiF^^(®) + In (l + e2y^^(®) + egF^^^®) + e.Y'"'^''^) 

where S{Q) is given by: 

5'(@) . y' = y[( ^^^ ) yC/^n./^ia) + ( ^12 ] y(?»2i,/i22) _^ | ^13 \ Y{h-iuh32h 

l C21 / I C22 / I C23 / 



Using m) of ( [5.5| ) we have Bi = B4 = H2 — Hi and looking at the array A, 
we see that (i?i, ^4; Cj) = if and only if j 7^ 1. Thus we have 

a) {H2 — Hi,Cj) = if and only if j 7^ 1. Similarly, 

b) {B2;Ci)=0 implies {3H2-2Hi-H3;Ci) = (} 

c) C3 = mC2. FoT,{H2 - i^i; C2, C-s) = implies i/i = H2, 

if C2, C3 are hnearly independent, 
since, (the order of S{Q) is two.) 

Thus S{Q) is reduced, by these conditions to a 9 parameter system. 

1(9) = ei(e)r^^(®) + \n(i + J2 efc(e)y^^(®)) 

\ k=2 J 



and from (|575| ), ( |5.6| ) we can find i?fc(9), efc(9)(A; = 1, ... ,4). In the special 
case when 5'(9) comes from a second order differential equation, we have the 
additional conditions: 

5(6) : v' = yl(l) y'-"^ + ( cl ) ^""''" + ( c! ] y('"'"+2'»»>l 

writing S{Q) as a second order differential equation we get: 
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y" = C22y^'' {y'f + c^2y'''''^^{yT''^^ 5.7 

/(e) = eiy('^2i+i) + In (l + e2{yT~''''^ + e^y^^''''+^^ + e^y^'^''^^'^) 
where 

^ ^ (fe32-l)c22 „ ^ (fe32-l)^(ci2) 

'^l /121+1 ' 2 2(/l21+l)2c32 



_ (fe32-l)^ci2 p ^ (l-fe32)c22 

Setting /i2i = 0, h^2 = — 1, C22 = 2, C32 = — 3 in (|5.7|) and its first 



integral /(6), we get (|5.1| ) and its first integral /. 

The differential equation ( |5.7| ) and its first integral /(6), is, actually, a 
special case of the following: 

y" = c22y''''{y'f + c2zy''{y'f 5.8 

where a = {q - l)/i2i + g - 2, /? = /i32 + 1. 

For any integer g > 3, (|5.8| ) has a first integral: 

-fc+i){/i3i+i) 1 5.9 

fc=3 



J = ei2/('^^^+^) + In 1 + e2y'^^-^^^^ + ^ Cfcy^" 



Setting g = 4 in (p.8D, (p.9|), one gets ( p. 71) and its first integral !{&)■ 



6 The 2? X g Array (the algebraic case). 

In the following we describe the role that the general p x q array plays in 
determing the necessary and sufficient conditions, that the system S given 
by (|1.3 ), must satisfy, for the system to have an algebraic first integral given 
byO- 

Let / = J2t=i ^kY^'" be a first integral of the multinomial system given 
by ( p. .31 ). Making use of the derivative formula i) of ( p.l| ) we get: 

= ELi(E,^=ie.(S.;C,)F^'=+^0 = 
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In the set of g^ vectors Bj. + Hj, leave out all vectors such that (5^; Cj) = 
0. Let the remaining vectors form a set of p distinct vectors Ei, . . . , Ep, then: 



/' = E[ E e,(i?,;C„J]y^'=0 



6.1 



-1 Bk+Ha-—Ei 



Since i?, (i = l,...,p) are distinct, equation ( |6.1| ) yields a system of]? 
homogeneous linear equations: 



fc=i 



6.2 



which ei, . . . , Cq must satisfy. The p x g array A is a pictorial representation 
of (B.ll) and is defined as follows: 



A = (A 



ik) 



1, . . . ,p;k = 1, . . . ,q), where: 



■^ik 



H^ if Bk + H^ = Ei and (5^; CJ ^ 
if no such Ha exists 



6.3 



The symbol stands for the empty spot. 
Definition(|6.3|) implies: 



i) Ha appears in the A;*^ column of A if and only if 
{Bk, Ca) + 



ii) if Ha,Hf3 appears in the same row in columns j,k, 
respectively, then: 

Bj + Ha = Bk + Hi3 



} 6.4 



We call columns: j, k linked if ii) of ( |6.4|) is satisfied and we set 



Ljk — Ha — Hp 
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We call j, k oi A connected if there exist columns. Pi, . . . , Pg such that 
j = Pi, k = Pa and columns P^, Pa+i (« = 1, • • • , s — 1) are linked, and we 

set 

s-l 
a=l 

e.g. In the array: 

[ H^ H, j 

columns 1,2 are linked and L12 = H^ — H^; columns 2,3 are linked and L23 = 
Hyj — H^; columns 1,3 are connected and L13 = Lu + L23 = Hu + H^ — 2Hy. 
Note, that when the array A is connected, Ljk is defined and 

Ljk = Bk — Bj for all 1 < j < /c < g 

This implies that the difference of any 2 exponent vectors of / is a linear 
combination of the exponent vectors of the system S. 

It is not difficult to show that when A is not connected (i.e. when there 
exist at least two column of A which can not be connected) then A is the 
array of an integral I oi S such that I = I1 + I2 where Ji and I2 are both first 
integrals of S, [1] . From now on we shall assume that the array A is connected. 
It follows from the definition of A that interchanging columns j, A; of A is 
equivalent to interchanging the exponents Bj, Bk of / and interchanging rows 
i,j of A is equivalent to interchanging Ei,Ej of (|6.1|) . Thus we identify all 
arrays that can be obtained from one another by an interchange of rows and 
or columns. 

Along with A we define the p x q matrix 

MA = (ttifc {i = l,...,p;k = l,...,q) 
where 



O-ik — {Bk] Co) if Aik — Ha 

ttik = if Afc = 



6.5 
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There are two kinds of arrays normal and abnormal. An array A is called 
normal if for every H^ that appears in A, there exist at least one column of 
A which does not contain Ha- An array A is called abnormal if there exists 
at least one H^ which appears in every column of A. It is remarkable that 
when A is normal we can compute [Bk] Ca) without knowing what the B^s 
are. For, let Ha fail to appear in column j of A then, for any k ^ j we have 

(Bk] Ca) = (Bk — Bj] Ca) + (Bj] Ca) 



Thus, M(A), which is the matrix of the system of homogeneous linear equa- 
tions given by ( |6.2| ), can be computed, when A is normal. 

Let A be an array of the exponent vectors of a system S, given by ( [L.3|) . 
The following is a set of necessary and sufficient conditions that A must 
sotisfy for A to be an integral array (i.e. there exists a first integral I oi S 
such that A satisfies (|6.tJ|)). 
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a) Each row of A must contain at least two distinct H^s ' 
of S and no H^ may appear more than once in any row or 
column of A. 

b) The Ljk (1 < i < A; < g) are well defined and do 
not equal (0 = the zero vector). 

c) If Ha fails to appear in columns j, k of A, then {Ljj^; Cq = 

d)If Ha appears in column j but fails to appear in col- 
umn k, then {Ljk] Co) 7^ 0. 



e) Rank of any q — 1 columns of M(A) equals rank of 
M(A) =g-l. 

f) Bi satisfies the linear system: 

{Bi, Ca) = if Ha fails to appear in column 1. 

{Bi, Ca) = {Lji, Ca) if Ha appears in column 1, but fails 

to appear in column j (the exis- 
tence of such a "/' is guaranteed 
by the fact that A is normal). 



> 6.7 



Condition a) is necessary, since /' = implies that the coefficients of 
Y^^ (i = 1, . . . ,p) must vanish, thus if only one Ha appears in row i then 
the coefficient of Y^^ is = contrary to assumption. Also, 

Ha = Aij = Aik implies Bj = B^ and 

Ha = Ak = A,jk implies E., = Ej 

condition b) is a strong restriction on S e.g. if 

I Hi H2 \ 



A 



H-t Ha 



then L12 = Hi — H2 = H3 — H4 yields a linear relation among the exponents. 
Also, Ljfc = implies Bj = B^. 
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Conditions c),d) are implied by Lj^ = B^. — Bj, thus if H^ does not 
appear in either column, we have {Lj^, Co) = {B^; Co) — {Bj, Co) = + = 
0. Similarly if Ha appears in column k but not in column j, [Ljk] Co) = 
{Bk] Co) — {Bj, Co) = {Bk, Co) — 7^ 0. 

Condition e) is the requirement that the system of homogeneous linear 
equations in e^ {k = 1, . . . ,q), given by 

^aifcefc = i = {l,...,p), 
fc=i 

has a solution such that Y[T=i ^fc 7^ 0. 

Condition /) follows from (|6.3| ) and (|6.6|) by setting A; = 1. 

We now show that conditions a), . . . , f) are also sufficient for the con- 
struction of a first integral / = I]fc=i CfeF^* of 5" given by (|1.3| ). 

Let A = [Aik), M(A) = (uik) (z = 1, . . . ,p ; A; = 1, . . . ,g), where atk = 
if Aik = 0, o,ik = {Ljk;Ca) if Aik = Ha and j is such that Ha does not 
appear in column j of A. Let A, M(A) satisfy conditions a), . . . , /) then 
we constant a first integral / as follows: 

i) Let Bi be a solution of the linear system as in f). 

ii) Set Bk = Bi + Lik 

in) Set Ci, . . . , Cq to be a solution of the homogeneous linear 

system: 

J2k=i (^ik^k = 0, (a solution is guaranteed by e). 

iv)SetI = El=iekY''>'. 

We first show that {Bk, C„) = if and only if Ha does not appear in the 
k^'^ column of A. By i) this holds for A; = 1, let fc > 1. Let Ha not appear in 
the k^^ column of A and let Ha also fail to appear in the first column of A, 
then 

(Bk, Ca) = {B,; Ca) + {L.k, Ca) =0 + = 

{Bi, Ca) = follows from the definition of Bi given by i), {Lik, Ca) = 
follows from condition c) that A satisfies. Let Ha fail to appear in column k 
but appear in column 1, then {Bi, Ca) = {Lk, Ca) by i), so that 

(Bk] Ca) = {Bi, Ca) + [Llk] Ca) = [Lkl] Ca) + (-^Ifc; Ca) = 0. 

Now, let Ha appear in column k but fail to appear in column 1, then 

{Bk, Ca) = {Bi, Ca) + {Lik, Ca) = + {Lik, Ca) 7^ 0, 
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by condition d) that A satisfies. 

Finally, if H^ appears in column k and column 1, but fails to appear in 
some column j (such a j exists because A is normal) then 

(-Sfe! Co) = (Bi, Co) + [Lik] Co) = {Lji, Co) + {Lik] Co) = [Ljk] Co) 7^ , 

by condition d). 

This proves our assertion for all 1 < A; < q'. Let Ha = A^, Hp = Ai^, we 
show that 

Bj + Ha = Bk + Hfs = Ei (^ = l,...,p) 

For, 

Bk-B, = B^ + L^k-{Bi + L^,) {hyii)) 
= Lik + Lji = Ljk 

and Ha — Hp = Ljj. (by definition of Lj^ when columns j, k are linked) Thus 

Ha — Hr = B}^ — Bj and 

Ha + Bj = Hj3 + Bk = Ei 

We can, now, show that 

fc=i 
is a first integral of S. For, 

fc=l k=l jr = l 

But {Bk, Cj) = for any Hj which does not appear in column k, therefor 
the sum of the {Bk, Cj)'s is restricted to summing along the columns of A. 
If instead we sum along the rows of A, we get: 

k Q 

4=1 k=l 

For each Hj in A, let aj be the column of A in which Hj fails to appear, 
then: 
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{Bk] Cj) — {Bk — Baj] Cj) + {Baj, Cj 

= {Laj,k; Cj) + = {Laj,k; Cj) as in (|6 
= aik by (|6 

Thus 

I' = T.iT.^ka^k)Y'^^=0 

i=l k=l 

by the choice of e^ (A; = 1, . . . , g) in in). 
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